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Abstract. We present a new grid of LTE model atmospheres for weakly magnetic (B 10 10 G) neutron stars, using 
opacity and equation of state data from the OPAL project and employing a fully frequency- and angle-dependent 
radiation transfer. We discuss the differences from earlier models, including a comparison with a detailed NLTE 
calculation. We suggest heating of the outer layers of the neutron star atmosphere as an explanation for the 
featureless X-ray spectra of RX J1856.5-3754 and RX J0720.4-3125 recently observed with Chandra and XMM. 
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1. Introduction 

Modern X-ray observatories detect the thermal surface 
emission of a number of neutron stars. In analogy to the 
classic model atmosphere analysis of normal stars, such 
observations permit the direct measurement of fundamen- 
tal neutron star properties, such as their effective tem- 
peratures, atmospheric abundances and surface gravities. 
The latter point is especially of great importance, as an 
accurate measurement of the surface gravity of a neutron 
star is (with the distance known) equivalent to a measure- 
ment of its mass/radius ratio. The observational confirma- 
tion/rejection of the predicted neutron star mass-radius 
relations is a fundamental test of our understanding of 
the physics of matter above nuclear densities. 

The first neutron star model at mospheres invo lving re- 
alistic opacities were computed by Romani (1987 , hence- 
forth R87), using atomic data from the Los Alamos 
Opacity Library and employing a simple angle-averaged 



tion transfer as R87. Contemporaneously, a similar set of 
low-field atmospheres, partially based on the same atomic 
OPAL data, but employin g a more sophisticat ed radiation 
transfer, was presented by Zavlin et al. (1996 , henceforth 
ZPS96). The RR96 and ZPS96 models, broadly confirm- 
ing the results of R87, were applied to X-ray observations 
of presumably low-field neutron stars, i.e. millisecond pul- 
s ars (RR96, Zavlin fe Pavlov 1998 ), isolated neutron stars 
(Pavlov et al. 1996), and transiently accreting neutron 



stars in quiescent LMXBs (e.g., Rutledge et al. 1999 , 2000 
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a inajui result; itsr cur 
the thdiinal uinissiuii uf iieuLiun sLais diffeis substantially 

from a Planck spectrum. For low-metallicity (helium) at- 
mospheres, the emitted spectrum is harder than the corre- 
sponding blackbody spectrum. The spectra emitted from 
high-metallicity atmospheres (carbon, oxygen or iron) are 
closer to a blackbody distribution, but show strong ab- 
sorption structures in the energy range observable with 
X-ray telescopes. 

Rajagopal & Romani (1996, henceforth RR96) com- 
puted hydrogen, solar abundance, and iron model atmo- 
spheres for low-field neutron stars, using improved opac- 
i ty and equation of sta te data from the OPAL project 
([glesias fe Rogers 1996), but employing the same radia- 



In the case of strongly magnetic neutron stars (B ~ 
10 12 — 10 14 G), both the observation and the computation 
of the thermal surface emission is significantly more diffi- 
cult, as non-thermal magnetospheric radiation dilutes the 
surface emission and as little atomic data is readily avail- 
able (for a review o f the pro perties of matter in strong 
magnetic fields, see Lai 2001 ). Magnetic models for hy- 
drogen and iron atmospheres were computed by Pavlov 



et al. (1995) and Rajagopal et al. (1997), respectively. 
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The wealth of high-quality neutron star observations 
expected from Chandra and XMM justifies the compu- 
tation of modern model atmospheres. In this paper, we 
present new model atmosphere grids for low-field neutron 
stars that overcome a number of shortcomings and errors 
in the RR96 and ZPS96 calculations and that are made 
available to the community. 

2. The model atmospheres 

The new grid of neutron star atmospheres was computed 
combining a fully f requency- and angle-d ependent radia- 
tion transfer code ( Gansicke et al. 1995 ) and the OPAL 
opacities from RR96. The employed radiation transfer 
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is part of a model atmosphere code that has previously 
been used for the computation of white dwarf atmospheres 
which have been widely applied to observations of cata- 
clysmic variables and white dwarf/main-scquence binary 
stars. 

The model atmospheres were created under the follow- 
ing assumptions. (1) Plane-parallel geometry. The atmo- 
sphere of a neutron star has a vertical extension of a few 
centimeters at most, compared to a radius of ~ 10 km. Its 
curvature is hence completely negligible. (2) Hydrostatic 
equilibrium. For the extreme surface gravity of neutron 
stars and the relatively low temperatures that we con- 
sider here, the atmosphere is static. The situation is differ- 
ent in accreting neutron stars during X-ray bursts, where 
the nuclear burning atmosphere considerably expands on 
short time scales. (3) Radiative equilibrium. The atmo- 
sphere contains no source of energy, but acts only as a 
blanket, through which the thermal energy of the under- 
lying core leaks out. In Sect.||, we will discuss how aban- 
doning this assumption will impact the spectrum of the 
neutron star. (4) Local thermal equilibrium (LTE). For 
the high densities and the low temperatures in the consid- 
ered atmospheres collisional ion-ion interactions dominate 
over interactions between matter and the radiation field 
throughout most of the atmosphere. Deviations from LTE 
in the outermost tenuous layers of the atmosphere may 
somewhat affect the depth of the absorption lines, but 
considering the overall level of uncertainty in the atomic 
physics involved, and the relatively poor quality of the 
available spectroscopy, the assumption of LTE appears to 
be a reasonable approximation. Below, we discuss a quan- 
titative comparison to a NLTE model. 

With these approximations, the computation of a 
model atmosphere can be separated into two independent 
parts: calculating the structure of the atmosphere from 
the integration of the hydrostatic e qua tion (Sect. 2A) and 
solving the radiation transfer (Sect. |2~2| ). The temperature 
struct ure of the atmosphere. T(z). is a free parameter in 
the ov erall process and is adjusted iterativclv until radia- 
tive and hydrostatic equilibrium are satisfied (Sect. |2~3| ). 

2.1. Atmosphere structure 

The vertical scale of the atmosphere is transformed from 
the geometrical depth z to the optical depth r by dz = 
dr/px, with \ the total mean opacity (see below) and p the 
density. The structure of the atmosphere is then obtained 
by integrating the equation of hydrostatic equilibrium 



mosphere are isothermal and that the degree of ionization 
is constant. With \ = const, 



X(Tmin) 



(2) 



1 



X being itself a function of (T, P g ), we use x(r m i n 
as an estimate and solve eq. (^) with a Newton-Raphson 
method. From P g (r — T m i n ) we proceed with the integra- 
tion of eq. (0) to T max . 

This integration involves at each depth the evaluation 
of the equation of state (EOS) as well as the computa- 
tion of the absorption and scattering coefficients. For the 
EOS and the radiative opacities, we use the OPAL data 



(Iglesias & Rogers 1996) from RR96, that cover three dif- 
ferent chemical compositions: hydrogen, solar abundances, 
and iron. The opacities are tabulated as a function of the 
temperature T, R — p/T§ with Tq the temperature in 
10 6 K, and u — E^/ksT, with E 7 the photon energy (for 
additional details on the OPAL opacity tables, see RR96). 
For a given depth in the atmosphere, defining T and p 
and at a given energy, the radiative absorption coefficient 
K ra d (T, p, E) is calculated from the opacity tables using a 
bilinear interpolation in log T and log R and a linear in- 
terpolation in E. In the construction of the model atmo- 
spheres, the radiative absorption coefficient K ra d(T, p, SE) 
for a finite energy interval SE is given by the harmonic 
mean of n individual evaluations of K ra d(T, p, i?„), with 
E n within 5E, and n sufficiently large to sample well the 
resolution of the OPAL tables. Thomson scattering is an 
important source of opacity only in low-Z atmospheres 
at high temperatures and high energies. The OPAL tables 
implicitly include Thomson scattering as radiative opacity. 
K ra d is, hence, the sum of true absorption plus Thomson 
scattering. We will come back to this issue in Sect. |2.2| . 
Finally, we compute a radiative Rosseland mean opacity 
kross from the energy dependent K la d(E). 

Following RR96, we include a conductive opacity 



K con = (2 5 x 10 4 /n e ) x (Z a /A)(T 7 1/a /p)cm a g- 1 QCox k 
Giuli 1968) ) to account for energy transport by electron 
conduction, with n e the number of electrons per a.m.u., 
A the average atomic mass, Z = An e the average ionic 
charge, and T7 the temperature in 10 7 K. The radiative 
(Rosseland) opacity and the conductive opacity are har- 
monically added. The total opacity is, hence, given by 

X = ( K Ross + (%on) • 

Once eq. (Q) is integrated to r max , a complete macro- 
scopic description of the atmosphere structure is at hand. 

2.2. Radiation transfer 

For the assumption given above, the radiation transfer 
equation takes the form 



dPg 
dr 



9_ 

X 



(1) 



where P g is the gas pressure, from an outer boundary r min 
to an inner boundary r max . Considering the extremely 
high surface gravity encountered in a neutron star atmo- 
sphere, we neglect the effects of radiation pressure. We 
compute a starting value for the integration of eq. (|l|), 
P g (r — T m ; n ), assuming that the outer layers of the at- 



d ly(p) 
dT v 



I v {p) - S v , 



(3) 



with p, = cos 9 the cosine of the angle between the normal 
to the atmosphere and the direction of the considered ra- 
diation beam, t„ the frequency-dependent optical depth, 
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Fig. 1. Emergent model spectra (unredshifted Eddington flux) for M ns = 1.4-Mq, R lls = 10km neutron stars. 
From left to right: hydrogen, solar abundance, and iron atmospheres. From top to bottom in each panel: logT ns = 
6.5, 6.3, 6.1, 5.9, 5.7, 5.5, 5.3. Plotted as grey dashed lines are blackbody spectra for the corresponding temperatures. 



I„{p) the frequency- and angle-dependent specific inten- 
sity, and S v the source function. 



As mentioned in Sect. 2.1, Thomson scattering is im- 
plicitely included in the OPAL opacity tables. Considering 
that the contribution of Thomson scattering to the total 
opacity is negligible in high-Z atmospheres, we restrict an 
explicit treatment of Thomson scattering to our hydro- 
gen models (Sect.pl). We compute Thomson scattering, 
a e = N e 6.65 x 10~ / p, but taking into account the cross- 
section reduction by collective effects (Boercker 1987), 
with N e the electron density in cm' 3 , calculated under 
the assumption of full ionisazion, and the true frequency- 
dependent absorption «„ = K ta d(y) — n e , where K ra d(^) 
is the radiative frequency-dependent OPAL opacity from 



Sect. 2.2 



With the isotropic and, to first order approximation, 
coherent Thomson scattering term, the source function is 
given by 



From the resulting angle-dependent specific intensities, 
the flux at the surface of the atmosphere is computed, 



F v = pl v (p)dp 



Radiative and conductive equilibrium requires 



F v du + F con = 



(5) 



(6) 



at each depth in the atmosphere, where a is the Stefan- 
Boltzmann constant, F con d = (16oT 3 ) / (3pK con )dT j dz is 
the conductive energy flux, as detailed in RR96, and T ns is 
the effective temperature of the neutron star. If the devia- 
tions from the equilibrium (|^) are larger than a given error, 
correction terms AT(t) are computed fro m linearization 
of the angle- integrated form of cq. (|^) (e.g. Feautrier 1964; 
Aucr fc Mihalas 1969|, |l970j; |Gustafsson 1971 



S v = 



-B v 



- J,' 



(4) 



with B v the Planck function and J v the mean intensity. 
We use as boundary conditions for eq. I(p < 0) = at 
T m j n and the diffusion approximation at r max . Equation (||) 
is solved using the well-documented Rybicki method (e.g. 
Mihalas 1978b - 

For the high-Z composition atmospheres (solar abun- 
dances and iron), we treat Thomson scattering as true 
absorption, directly using the OPAL K ra d- 



2.3. Model atmosphere computation 

The actual computation of a model atmosphere and the 
emergent spectrum is done in two separate stages. 

(i) Starting from the temperature structure of a grey 
atmosphere (|Chandrasekhar 1944 ), we iteratively con- 
struct (Sect. 2.1 and Sect. 2.2) an atmosphere structure 
that satisfies radiative and hydrostatic equilibrium. At 
this stage, we use an energy grid of 1000 logarithmic 
equidistant bins covering the range 10~ 3 4 — lOkeV. For 
each energy bin, 30 evaluations of the OPAL opacity tables 
are harmonically added for the computation of the radia- 
tive absorption coefficient. The atmosphere structure is 
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Fig. 2. Ultraviolet model spectra of non- magnetic neu- 
tron stars in the wavelength range observable with 
HST/STIS, redshifted for z = 0.306. Solid lines: iron com- 
position. Dashed lines: hydrogen composition. From top to 
bottom: logT ns [if] = 5.7,5.5,5.3,5.1. The absolute fluxes 
are given for R ns = 10 km, d = 100 pc, and E(B — V) = 0. 



constructed on an optical depth grid with 100 points cov- 
ering tr = 10~ 6 — 5 x 10 3 (except for the coldest models, 
which were computed on a grid covering tr = 10~ 6 — 10 3 
due to limitations in the available OPAL data). The angu- 
lar dependency of the radiation field is sampled by three 
Gaussian points in fi. Radiative equilibrium (eq.|] is satis- 
fied to better then 10~ 5 at each optical depth) is reached 
after ^5 — 15 iterations. Generally, the hydrogen atmo- 
spheres converged the slowest, as the steep drop of the 
opacity oc u makes the atmosphere transparent for high- 
energy photons, thus radiatively coupling the deep hot 
layers with the surface layers. Electron conduction carries 
a few percent of the total flux at the largest optical depths 
(Rosseland optical depths tr oss > 100) in the coldest iron 
and solar atmospheres. For log(T ns ) > 5.75, the conduc- 
tive flux is less than 1%, which is in good agreement with 
the estimates of ZPS96. Including the conductive opacity 
in our atmosphere models has, hence, no noticeable effect 
on the emergent spectra. 

(ii) From a given atmosphere structure, we recom- 
pute the emergent spectrum solving the radiation trans- 
fer (Sect[2.2|) on a much finer energy grid (10000 bins, 10 
OPAL evaluations per bin) , and using 6 Gaussian points in 
/i to resolve the angular dependence of the specific inten- 
sity. The higher number of angular points has practically 
no influence on the angle-averaged flux from the neutron 
star surface, but may be used to account for limb darken- 
ing in the computation of spectra from neutron stars with 
non-homogeneous temperature distributions. 



2.4. The model spectra 

We computed grids of model atmospheres for the three 
different compositions covered by the OPAL opacity/EOS 



tables of RR96: hydrogen, solar abundances, and iron. We 
restricted ourselves to the canonical neutron star configu- 
ration, Af ns = 1.4 Mq and R ns = 10 km, corresponding to 
a surface gravitational acceleration of log g = 14.386. The 
three grids cover the temperature range logT ns = 5.1 — 6.5 
in steps of 0.05, i.e. a total of 29 spectra per abundance 
grid.[] 

Figure [l] shows the emergent spectra for the three 
different compositions. As already evident in the previ- 
ous neutron star atmosphere calculations (R87, RR96, 
ZPS96), the model spectra differ significantly from black- 
body distributions. The pure hydrogen models show a 
strong flux excess over the blackbody distributions at en- 
ergies above the peak flux. The free-free and bound-free 
opacity in these atmospheres drops off rapidly at high en- 
ergies, leaving Thomson scattering as the dominant inter- 
action between the atmosphere matter and the radiation 
field. As a consequence, the atmosphere is highly trans- 
parent to the hard X-ray photons from deep hot layers. In 
contrast to the hydrogen models, the iron and solar abun- 
dance models are overall closer to the blackbody distri- 
butions because of the milder energy dependence of their 
opacities. They show, however, a substantial amount of 
absorption structures, lines and edges, which are especially 
pronounced in the ~ 0.1 — lOkeV range well observable 
with most of the present and past X-ray satellites. 

Noticeable absorption structures with equivalent 
widths of ^ A are present in the iron spectra of moder- 
ately cool neutron stars also in the ultraviolet. The (non- 
magnetic) hydrogen models contain only the extremely 
pressure broadened Lya line. In a magnetic field Lya 
splits in three Zeeman-components, with the er + /~ com- 
poncnts asymmetric ally shifted by ~ 100 A for B w 10 9 G 
( Ruder et al. 1994 ). Because the hydrogen in these hot 
atmospheres is ionised to a large extent, the Lya Zeeman 
components are expected to be rather weak, but they may 
possibly be detected with future large aperture ultraviolet 
telecopes, allowing a direct measurement of the magnetic 
field strength. Figure|| compares the redshiftedQ iron and 
hydrogen spectra in the far ultraviolet. 

Figure H shows the angle-dependent intensity 
for a T ns = 10 6 K iron atmosphere. It is apparent 
that the emission from the neutron star surface is 
highly anisotropic, as previously discussed by ZPS96. 
Calculations of the emission of a neutron star with a non- 
homogeneous surface temperature distribution must take 
this anisotropy into consideration. An application of the 
angle-dependent intensities will be discussed elsewhere. 



1 The angle-averaged fluxes are available from CDS, for the 
angle-dependent specific intensities, or models covering addi- 
tional temperatures/surface gravities, please contact the au- 
thors. 

2 2 = (1 - 2GM ns /R ns c 2 )- 1/2 - 1 = 0.306 for an assumed 
Rns = 10 km and M na = 1.4 M©. 
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Fig. 3. The angle-dependent specific intensity of a 

T ns = 10 6 K iron model for p = 1.0, 0.7, 0.4, and 0.1. 

2.5. Comparison to earlier neutron star atmosphere 
models 



2.5.1. |Rajagopal & Romani (19961 ) 



The new model spectra agree quite well with the mod- 
els of RR96. The only systematic differences are found 
for the hydrogen models, where the new models have a 
significantly lower flux at energies il 1 keV. This differ- 
ence is due to the inclusion of Thomson scattering in the 
new models, which increases the opacities at large optical 
depths in the hydrogen atmospheres. For the heavy ele- 
ment atmospheres, this effect is negligible, as the opacity 
is dominated by bound-bound and bound-free absorption. 

Small (5-15%) differences in the emergent flux of 
the heavy element atmospheres are found in the low- 
temperature (logT ns < 5.75) models in the energy range 
1-10 keV. The specific intensities at these high energies 
are 15 — 25 orders of magnitude below the peak intensi- 
ties, and we believe that the differences between the new 
models and those of RR96 are due to the different numeri- 
cal treatment of the radiation transfer. At higher energies, 
the agreement between the two different model genera- 
tions is better than 1%, except in strong lines where dif- 
ferent energy sampling naturally leads to somewhat larger 
discrepancies. 



2.5.2. iZavlin et al. (1996ft 



ZPS96 computed grids of non-magnetic neutron star LTE 
model atmospheres for three different compositions, hy- 
drogen, helium, and iron. The hydrogen and helium opac- 
ities and non-ideal EOS employed in these models were 
calculated by ZPS96 using the occupation probability for- 



malism described by Hummer & Mihalas (1988). The iron 
models were computed using the same OPAL data used 
also by RR96 and in the present paper. The iron models 
show a significant difference both in the overall shape of 
the continuum as well as in the details of the absorption 



Fig. 4. Comparison of the T n; 



200 000 K pure iron 



NLTE model spectrum from Werner fc Deetjen (2000 ) 
(gray curve, unredshifted Eddington flux) with a corre- 
sponding LTE model spectrum from our new atmosphere 
grid (curve). 



features. We traced this to an erroneous use of the opac- 
ity grid in ZPS96. This has been confirmed (Pavlov, priv. 
comm.)^]. The ZPS96 hydrogen model is somewhat harder 
than both our new spectrum and that of RR96. The flux 
excess of the ZPS96 model (or the flux deficiency of our 
model) is not fully understood, but likely stems from small 
differences in the absorption/scattering opacities. 



2.5.3. [Werner & Deetjen (200q ) 

'Werner & Deetjen (2069) presented the first NLTE cal- 
culation for a neutron star atmosphere using bound-free 
opacities from the Opacity Project ( Seaton et al. 1994 ) 
and explicitly treated line blanketing by millions of lines 
from the Kurucz & Bell (1995) line list. Figure^ compares 
the T ns _ 200 000 K pure iron NLTE model from |Wcrnci 
|fc Deetjen (2000 ) with a corresponding model from our 
LTE calculations. While the NLTE model clearly shows 
more detail because of the much larger number of consid- 
ered transitions, the overall agreement between the two 



models is very good. Werner & Deetjen (2000) quote a 
flux difference between their LTE and NLTE models of 
~ 10 % in the line cores, and much less in the contin- 
uum. For higher temperatures, the NLTE effects might be 
stronger, but so far no detailed atomic data are available. 

3. Beyond "classic" stellar atmospheres 

Part of the ongoing effort in computing neutron star atmo- 
spheres has been carried out for the an alysis of X-ray ob - 
servations of isolated neutron stars (e.g. Pons et al. 2002). 



3 After we alerted G. Pavlov to their errors in employing the 
opacities, their group has computed revised models in reason- 
able agreem ent with other work. A n example revised spectrum 
appears in (Pavlov & Zavlin 200C). 
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Fig. 5. Temperature structures T(t) in a locally heated 
neutron star atmosphere. An energy amount of 0.035crT ns 
(0.82a-r n s 4 ) is deposited at an optical depth Td op = 0.01 
(Tdep = 1). This energy injection is spread by a Gaussian 
distribution with a width of logr = 0.3, as illustrated in 
the figure, avoiding unrealistic (and numerically problem- 
atic) step functions in the T(r) structure. The undisturbed 
temperature run is plotted as solid line, the atmosphere 
heated at r = 0.01 (r = 1) as dashed (dot-dashed) line. 



However, recent high-resolution X-ray spectroscopy of the 
brightest isolated neutron star candidates RX J 1856. 5- 



3754 (Burwitz et al. 2001) apparently excludes all the neu- 



tron star atmosphere models presented so far: while the 
overall spectral energy distribution (SED) rules out \ow-Z 
atmospheres, the featureless XMM and Chandra spectra 
are inconsistent with the strong structures expected from 
a heavy element atmosphere. 

As in all previous calculations (R87, RR96, ZPS96, 
Werner & Deetjen 2000), the model atmospheres pre- 
sented here assume radiative equilibrium. It is worth re- 
membering that when the model fits suggest that the flux 
is dominated by hot polar caps with a small fraction of the 
full neutron star surface area, some sort of local heating 
is likely implicated. For instance, anisotropic interior con- 
ductivities can produce smooth variations in the surface 



effective temperture; e.g. Heyl & Hcrnquist (1998). For 
slow, long P neutr on stars, local heating can b e caused by 
low M_ accretion. Zampieri et al. (1995 ) and Zane et al 



(2000^ have produced emergent spectra for ionized H stars 
accreting at low rates. With proton stopping depths of 
~ 20 g cm -2 , much of the energy is deposited fairly deep 
in the atmosphere, but shocks in the accretion flow can 
apparently heat the outer layers, providing excess flux on 
the Rayleigh- Jeans tail. For an active magnetosphere, lo- 
cal heating without accretion arises from precipitating e^- 
(Arons 1981) or from illumination by energetic (polar cap 
or outer gap) photons. 

An analogous situation - pole caps heated by accre- 
tion - is well-documented in polars, a subtype of cata- 
clysmic variables in which a strongly magnetic white dwarf 




-0.5 
log 10 (E [keV]) 

Fig. 6. Emergent spectra from a heated T ns = 10 6 K iron 
atmosphere. Top panel: 0.82crT ns 4 are deposited at Td cp = 
1, equivalent S = 3.2 x 10~ 3 gcm -2 . The spectrum of 
the heated atmosphere and of the undisturbed atmosphere 
are plotted as black and as gray lines, respectively. The 
dashed line shows a blackbody spectrum corresponding 
to the effective temperature of the heated atmosphere. 
Bottom panel: 0.035crT ns 4 are deposited at deposited at 
Tdcp = 0.01, equivalent £ = 2.2x 10~ 4 g cm~ 2 . 



accretes from a late type donor star. In theses systems, the 
white dwarf atmosphere near the magnetic poles is heated 
by strong irradiation with cyclotron radiation and thermal 
bremsstrahlung from a stand-off shock. As a consequence 
of this accretion heating, the spectra of these heated pole 
caps are almost devoid of spectral structures. The Lyman 
lines, typically very strong absorption features in the ul- 
traviolet spectra of white dwarfs, are almost completely 

and also the X-ray 



flattened out (Gansicke et al. 1998) 
spectra of polars are very close to blackbody distributions, 
unlike the expected emission from a high-gravity photo- 
sphere strongly enriched with heavy elements (Mauchc 



1999). While the physical processes involved in heating the 



pole caps of accreting white dwarfs and accreting neutron 
stars differ markedly, the observed phenomenology is very 
similar: RX J0720. 4-3125, another bright isolated neutron 
star shows a quasi-sinusoidal X-ray light curve suggesting 
the presence of large heated pole caps, whereas the X- 
ray spectrum contains no significant absorption structure 
(Haberl et al. 1997; Paerels et al. 2001) reminiscent of 



the ultraviolet observations of the heated white dwarf in 



the polar AM Her ( Gansicke et al. 1998 ). It appears likely, 
therefore, that heating effects may affect the line spectra 
of heavy element neutron star atmospheres. 
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As a simple illustration, we ignored the details of par- 
ticular physical heating effects and computed a number 
of spectra under the assumption that an energy aaT ns 
is deposited at a characteristic depth in the atmosphere, 
T dcp (Fig. |5|). In practice, the emergent spectrum of a 
"heated" neutron star is computed by (1) modifying the 
thermal structure of an undisturbed atmosphere T(r) to 
account for the deposited energy, (2) recomputing the at- 
mosphere st ruct ure based on the modified temperature 
profile (Sect. ^Tl) , and finally solving t he r adiation tansfer 
within the "heated" atmosphere (Sect. 2/2). Figure^ shows 
the emergent fluxes from a T ns = 10 6 K iron atmosphere 
for two different sets of parameters, a — 0.82 and Td op = 1, 
equivalent to S = 3.2 x 10~ 3 g cm -2 , and a — 0.035 and 
T dcp = 0.01, equivalent to S = 2.2 x 10~ 4 g cm~ 2 . 

Two effects are evident. The flatter temperature in 
the outer layers of the heated atmospheres supresses the 
strong absorption edge around 1 keV, returning the Wien 
tail colour quite close to that of a simple blackbody. 
Further, the equivalent widths of the lines are, of course, 
strongly affected. Heating at t ~ 1 strongly supresses the 
< keV absorption features; surface heating in fact drives 
the line features on the Rayleigh- Jeans tail into emission. 

We stress that this illustrative approach is by no means 
a self-consistent model of a heated neutron star atmo- 
sphere, as we include neither any physical assumption on 
the actual heating mechanims, e.g., there are currently no 
plausbile physical processes known which are capable of 
heating such shallow layers, nor a proper energy balance 
of the heating/cooling processes. It is merely intended as 
a possible explanation for the apparent paradox that the 
SEDs of isolated neutron stars (in particular RX J1856.5- 
3754) are best-fitte d with hea vy-element model spectra 



( Pavlov et al. 1996 ; this paper; Burwitz et al. 2001| ; |Pons 
et al. 2 [302| ), whereas their X-ray spectra show at best very 
weak spectral structures. 

More detailed future work will need to carefully exam- 
ine the efficiency of the various possible heating mecha- 
nisms, in particular to establish the depth-dependent en- 
ergy deposition, and to account properly for the energy 
balance in the neutron star atmosphere. 

4. Conclusions 

We have presented a set of angle-dependant low field 
neutron star atmosphere spectra. These radiative and 
hydrostatic equilibrium atmospheres should be useful 
to researchers pursuing exploratory fits of soft X- 
ray/UV/optical data. However, it is important to inject 
a note of caution into the discussion. At a minimum, 
non-uniform surface temperature, combined with limb- 
darkening and gravitational focussing will have subtle, but 
important effects which can be modeled using the present 
grid. 

Many other physical effects can, of course, strongly 
perturb neutron star spectra. The effects on heavy element 
absorption lines and edges can be particularly strong. 
Magnetic fields characteristic of young radio pulsars will 



have a dramatic effect on the opacities and the emergent 
spectra. When considering the strength of the absorption 
lines in spectra from magnetic atmospheres with heavy el- 
ements ( Rajagopal et al. 1997| ) it is important to remem- 
ber that the line energies depend sensitively on the B-field. 
Even simple dipole B variation across a polar cap can shift 
line energies by ~ 10%, strongly decreasing the equiva- 
lent width of spectral features in phase-averaged spectra. 
Convection has also been considered (eg. RR96), although 
even small magnetic fields may supress convective trans- 
port and keep atmospheres close to the hydrostatic solu- 
tion. 

On the whole, while we feel that our atmosphere grid 
should be of use for fitting of observed data; complica- 
tions are certainly expected. In particular strong _B-ficld 
variations and surface re-heating can decrease the equiva- 
lent width of heavy element line features in neutron stars 
with active magnetosphers; this may explain in part the 
difficulty i n finding such features in early Chandr a/XMM 
data (e.g. |Paerels et al. 200l|; [Biirwitz et al. 200l|). 
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